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ABSTRACT 

Jets producing gamma-ray bursts (GRBs) are likely to carry a neutron component that drifts with respect to the 
proton component. The neutron-proton collisions strongly heat the jet and generate electron-positron pairs. We 
investigate radiation produced by this heating using a new numerical code. Our results confirm the recent claim 
that collisional heating generates the observed Band-type spectrum of GRBs. We extend the model to study 
the effects of magnetic fields on the emitted spectrum. We find that the spectrum peak remains near 1 MeV 
for the entire range of the magnetization parameter < £b < 2 that is explored in our simulations. The low- 
energy part of the spectrum softens with increasing eb, and a visible soft excess appears in the keV band. The 
high-energy part of the spectrum extends well above the GeV range and can contribute to the prompt emission 
observed by FenmfLPC[. Overall, the radiation spectrum created by the collisional mechanism appears to agree 
with observations, with no fine-tuning of parameters. 

Subject headings: gamma-ray burst: general — gamma rays: general — radiation mechanisms: non-thermal 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are produced by ultra- 
relativistic jets from short-lived and powerful energy sources, 
probably associated with black hole formation. Apart from 
the jet launching itself, a primary question of GRB theory 
concerns the emission mechanism of the burst: how does the 
jet emit the observed gamma-rays? 

The jet must develop an ultra-relativistic speed if its initial 
thermal energy is much larger than its baryonic rest-mass en- 
ergy. The thermal energy is dominated by radiation and can 
be released at the photospheric radius Ri, where the jet be- 
comes transparent. For sufficiently clean (baryon-poor) ther- 
mal jets, transparency occurs when radiation still carries most 
of the explosion energy. Such a "radiation-dominated" GRB 
has a high radiative efficiency and a spec trum that peaks near 
1 MeV (P aczynskilll986l:lGoodmanlll986). The shape of the 
spectrum in this case should be Planckian (lBeloborodovl20rTl 
hereafter B 11). 

The observed emission indeed peaks near 1 MeV in most 
bursts, however its spectrum is nonthermal, with an extended 
tail of high-energy emission. This has two implications: (1) 
the GRB photosphere is rarely radiation-dominated, and (2) 
some form of dissipation operates in the GRB jets and greatly 
broadens their radiation spectra. The free energy available 
for dissipation may be kinetic or magnetic. In particular, the 
relative motion of diff erent parts of the jet can be dissipated 
(e.g. lRees & MeszaroslfT99l iPaczynski & Xulll994 . It leads 
to collisionless shocks as well as collisional dissipation. 

Collisionless shocks are known as efficient accelerators of 
nonthermal particles in astrophysics. They offer a mecha- 
nism for generating nonthermal electrons in GRB jets, and 
their synchrotron emission may be associated with observed 
7-rays. The idea became popular, however it faces difficul- 
ties. First, recent ab in itio particle-in-cell simulations of col- 
lisionless shocks (e.g. 'S ironi & Spitkovskvl 120091 1201 Ih do 
not support this scenario. Their results suggest that nonther- 



mal electron acceleration by shocks is inefficient in GRB jets 
with expected transverse magnetic fields and magnetization 
parameter eb > 10^^- Secondly, the spectra of a large num- 
ber of GRBs observed by CG/?0/BATSE and FermilGBM in- 
struments are inconsistent with synchrotron emission, espe- 
ciall y in the fast-cooling regime that is expected in the model 
(e.g. iPreece et al.ll2000h . Some of the spectra show harder 
slopes than synchrotron emission can produce irrespective of 
the cooling regime. In addition, the synchrotron model does 
not predict the observed preferential position of the spectral 
peak, which is always in the MeV range for bright GRBs. A 
recent discu ssion of the internal -shock synchrotron model can 
be found in Daigne et al.l (1201 Ih . 

One can consider other phenomenological models of emis- 
sion from the optically thin region of the jet at r > i?^. 
In particular, one can consider synchrotron self-Compton 
radiation fro m continual ly and uniformly heate d plasma 
(IStern & Poutanen 20041 IVurm & PoutanenI l2009h . In this 
scenario, the optically thin plasma maintains a relativistic 
temperature, kT^ ^ rricC^, at which radiative cooling bal- 
ances the heating, and keeps radiating as the jet expands. 
While this model can produce hard spectra, it still does not ex- 
plain the preferential position of the observed spectrum peak. 
It also conflicts with the rece ntly observed high-energy spec- 
tra of GRBs Ce.g. lAbdo et ani2009bl IZhang et al.ll201 It) . 

A likely resolution of these problems is that the bulk of 
GRB luminosity is generated below the photosphere at r < 
rather than in the optically thin zone r > i?^. Then mul- 
tiple Compton scattering participates in the spectrum forma- 
tion and can naturally produce the observed spectra. Gen- 
erally, any subphotospheric heating leads to Comptoniza- 
tion of the (initially thermal) photons, and then radiation re- 
leased at the photosphere i?^ looks as nonth ermal while the 
peak of its spectrum remains near 1 MeV (Thompson 1994| 
Meszaros & Rees 2000b: Giannios 2006; Pe'eretal. 20O 



B eloborodo vl20 1 Ol hereafter BIO). A variety of Comptonized 
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spectra may be generated, depending on the assumptions of 
the model. Further progress can only be made when the phys- 
ical mechanism of energy dissipation is understood. It could 
be thermal or nonthermal, and its radial dependence is impor- 
tant for the emerging spectrum. 

Three heating mechanisms can operate in the subphoto- 
spheric region. Two of them (internal shocks and magnetic 
dissipation) rely of collisionless plasma processes, whose de- 
tails remain uncertain because of their complexity. The third 
mechanism (collisional heating) is straightforward and can be 
modeled from first principles. It was recently proposed that 
collisional heating shapes the GRB spectra (BIO). It will be 
the focus of the present paper 

Collisional heating is particularly str ong (and in- 
evitable) in jets that carry free neutrons (Derishev et alj 
[1999; Bahcall & Meszaros 2000; Fuller etal. 200rt 
Meszaros & Rees 2000a; Rossi et al. 2006; K oers & GianniosI 
12007). Before the complete decoupling of neutrons and pro- 
tons into two non-interacting components, the fading n-p 
collisions heat the protons to a mildly relativistic temperature. 
Importantly, the n-p collisions also create plasma that 
efficiently radiates its energy. The plasma is cooled by 
radiative losses and heated by Coulomb collisions with hot 
protons. The rates of all these processes are well known and 
their effect on the jet and radiation can be calculated without 
additional assumptions (BIO). In particular, one can show 
that the electron distribution function has two parts, thermal 
and nonthermal, and both Comptonize the radiation carried 
by the jet. 

Three features of collisional heating make it a promising 
emission mechanism: (1) It peaks at radii r ~ O.li?*, not too 
early and not too late to produce a bright photospheric emis- 
sion. (2) It has a high radiative efficiency, comparable to 50%. 
(3) The radiation spectrum produced by collisional heating 
can be accurately calculated from first principles. The striking 
result reported in BIO is that the spectrum has the Band-type 
shape with the high-energy photon index /3 ^ —2.5, which is 
consistent with observations. 

The formation of the radiation spectrum in a heated jet is a 
nonlinear problem that requires self-consistent calculation of 
the plasma and radiation behavior This problem was solved 
in BIO using a Monte-Carlo radiative transfer code in combi- 
nation with an iterative technique. Our present work has two 
goals: 

First, to provide an independent check of the result of 
BIO, we calculate the emitted spectrum using a different 
method: we solve the kinetic equations for the plasma and 
ra diation in the jet. We u se the numerical code described 
in IVurm & PoutanenI (120091 hereafter VP09). We develop a 
new version of the code to adapt it to GRB jets, calculate the 
emerging radiation spectrum, and compare the results with 
those of BIO. 

Secondly, the numerical models in BIO were limited to 
weakly magnetized flows, eb ^ 1, and neglected synchrotron 
emission. In the present paper, our models include syn- 
chrotron emission and self-absorption, and we explore the 
spectra produced by magnetized neutron-loaded jets for a 
range of < £b < 2. 

The paper is organized as follows. Section |2] describes the 
physical model and the setup of the simulations. Section|3]de- 
scribes our method of calculations. The results are presented 
in Sections |4]and|5]and discussed in Section|6] 



2. PHYSICAL MODEL AND SIMULATION SETUP 

The central engines of GRBs must be compact and hot 
objects, which are almost certainly neutron rich. The high 
density and temperature inevitably leads to / 3-equilibrium 
(Imsh enniketalJ II967I: iDerishev etall 1 1 9991; iBeloborodovl 
i2003i) . For all plausible parameters of the central engine, 
the equilibrium establishes a high fraction of free neutrons 
in the baryonic matter (this feature is seen, for example, in the 
accretion-disk model for GRBs, see Belobor odov 2008 for a 
review). Therefore, neutrons are generally expected in GRB 
jets. 

The presence of neutrons creates perfect conditions for 
strong collisional heating. In any variable baryonic jet, the 
inter-penetrating neutron and proton components inevitably 
develop relative motions (see BIO and references therein). 
The two components acquire different Lorentz factors Fn and 
F, and (rare) collisions between neutrons and protons dissi- 
pate enormous energy. The dissipation is strongest in those 
parts of the jet where F/Fn is highest. 

The variable jet consists of many causally disconnected 
shells. Different shells are accelerated to different Lorentz 
factors and emit different radiation, which is consistent with 
observed variability. In our simulations, we focus on the heat- 
ing history and radiation of one shell. The shell is launched 
from the central source and accelerated at the expense of its 
internal energy until it enters the coasting matter-dominated 
phase of expansion. At radius R^, the timescale for collisions 
between the neutron and proton components exceeds the ex- 
pansion timescale of the jet, and protons start to migrate rela- 
tive to neutrons, forming a compound flow with T T^- At 
this point, strong collisional heating begins and the radiation 
spectrum starts deviating from a blackbody. Our simulation 
starts at i?n and follows this evolution. We assume that F and 
F,i remain approximately constant till the end of the simula- 
tion (see BIO for discussion of this approximation). Note that 
migration of neutrons between different shells of a strongly 
variable jet can create compound flows with F » F„. 

i?n may be defined as the radius where the neutron flow 
becomes "optically thin" to collisions with protons. The cor- 
responding "optical depth" for this process at radius r is given 
by 



nncr„r 



Tn = 



Ln^n 



Rn 



47rTOnC'^rFi' 



(1) 



where is the comoving number density of the neutron flow, 
Ln = AnmnC^r^rf-^nn is its kinetic luminosity (isotropic 
equivalent) and t7n ^ 3 x 10^^^ cm^ is the effective cross- 
section for nuclear coUisions. Radius Rn is defined so that 

r„(i?„) = 1, 



Rn — 



LnCTn 



47rmnC'^F^ 



(2) 



We aim to calculate the evolution of electron/positron and 
photon distribution functions in the heated jet. The photon 
spectrum is controlled by electrons (and positrons) via Comp- 
ton scattering and synchrotron emission. Therefore, the key 
ingredient of any GRB model is how energy is injected into 
electrons. For the model studied in this paper, the electron 
heating is unambiguously determined by the collisional pro- 
cesses as described below (see BIO for detailed discussion). 

The electron heating comes from baryons, in two forms: 

I. Nuclear n-p collisions generate pairs whose energy 
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in the comoving frame of the jet is comparable to the 
pion rest-mass, -firicC^ ^ rriT^c^ « 140 MeV. The in- 
jection rate of these energetic particles is proportional 
to the rate of nuclear collisions, Ticoii — co'n'T-^nrrci, 
where Fid ~ ^(F/Fn + Fn/F) is the relative Lorentz 
factor of collisions. The injection rate is given by 
(BIO), 



1 Frcimp 

4 70^0 



3 r 



■ "coll 



F t 



(3) 



cxp 



Here n is the proton number density in the plasma rest 
frame and 



^cxp — ^ 



r 

cT 



(4) 



is the expansion timescale of the plasma, measured in 
its rest frame. The injection peaks near the Lorentz 
factor 7o w my^/mc w 300. We will approximate it by 
a Gaussian distribution centered at 70 — 300. These en- 
ergetic particles experience rapid radiative cooling and 
join the thermal population. 

2. The population receives energy from protons via 
Coulomb collisions. Note that protons have at least 
mildly relativistic random velocities — they are stirred 
by n-p collisions. The dominant majority of pairs 
are kept at a much lower temperature, because of 
Compton cooling. Coulomb collisions drain the en- 
ergy from protons to with the following rate (e.g. 
iGinzburg & Svrovatskii 1964) 



3 , , ctt moC ricn 
- m A 

2 Pp 



0.02., T^fL, (5) 

f-CXD 



where In A is the Coulomb logarithm, /3p ^ 1 is the 
random velocity of the protons stirred by nuclear colli- 
sions, ric = n_ is the number density of electrons 
and positrons, and 



TT 



(6) 



is a characteristic Thomson optical depth of the plasma, 
as seen by photons at radius r. The pair density ric 
is controlled by the process of injection (Equa- 
tion [3]), the subsequent pair-photon cascade in the ra- 
diation field, and annihilation, all of which must be 
self-consistently calculated. 

Equations (O and (|5]l determine the heating history of the 
expanding plasma shell and eventually the radiation spectrum 
emitted by this shell. The shell itself is a part of the (variable) 
jet. It has the following parameters: 

1 . Lorentz factor and kinetic luminosity (F and L for the 
plasma component and Fn and Ln for the neutron com- 
ponent). Note that the characteristic radius for a 
given shell is determined by its Fn and according to 
Equation (|2]l. 

2. Temperature of the thermal radiation at prior to the 
onset of collisional heating, TBB(^n)- The correspond- 
ing luminosity is Lbb — (4/3) caT^^T^AnRf-^ (here a 
is the radiation density constant). This parameter de- 
pends on L, F, and the radius at the base of the jet tq 
where our shell starts to accelerate. We fix tq = 10^ cm 



oc r ^ (assuming 
const, i.e. 



in our numerical models and calculate the correspond- 
ing TBB(^n) using the standard model of ad iabatically 
cooled outflow between rg and (see e.g. lPaczvnskil 
[1990; Piran et al. 1993). 

3. Magnetization £b- It is defined as the ratio of the mag- 
netic energy density Ub — B'^/Stt to the proper energy 
density of the plasma flow U = i/47rcr^F^ (which 
includes rest-mass energy). Both L^b and U are mea- 
sured in the comoving frame of the plasma. We assume 
that the magnetic field is advected from the central en- 
gine. Then the field must be transverse to the (radial) 
velocitjQ and scale with radius as B 
F « const). In this case eb = Ub/U 
does not change with radius. 

3. METHOD OF CALCULATIONS 

For our calculations we use a new version of the numeri- 
cal code developed by VP09. The code is designed to model 
the coupled evolution of the heated outflow and the radia- 
tion it carries. The evolution is tracked by solving the time- 
dependent kinetic equations for the particles and photons. 
The included interactions are Compton scattering, cyclo- 
synchrotron emission and absorption, photon-photon pair pro- 
duction and annihilation, and Coulomb collisions. We use the 
exact cross-sections or rates for all these processes. 

The original version of the code had one significant limi- 
tation: it used the "one-zone" or "leaking-box" approxima- 
tion. This approximation pictures a uniform plasma cloud of 
an optical depth tt with isotropic populations of particles and 
photons. It treats the loss of photons from the cloud using an 
escape probability instead of accurate calculations of radia- 
tion diffusion through the cloud (VP09). 

The leaking-box picture is not good for GRB jets. In con- 
trast to static sources, where radiation escapes the source on 
timescale ~ ttR/c, GRB radiation remains embedded in the 
relativistic jet at all radii of interest. It evolves with radius 
according to the radiative transfer equation, which is in seri- 
ous conflict with the one-zone approximation. For example, 
the true angular distribution of radiation in the jet comoving 
frame is far from being isotropic, even in the subphotospheric 
region (Bll). The isotropic approximation of the one-zone 
model becomes invalid when the optical depth tt decreases 
below ^ 10. Furthermore, the large free path of photons near 
the photosphere leads to mixing of radiation emitted by dif- 
ferent parts of the photospheric region with different Doppler 
shifts. This mixing has a significant effect on the local photon 
spectrum and the spectrum received by a distant observer 

For these reasons the simplified one-zone treatment has to 
be abandoned in favor of proper radiative transfer calcula- 
tions. The kinetic equation for photons in VP09 is replaced 
by the transfer equation as described below. 

3.L Radiative transfer 

The equation of radiative transfer in an ultra-relativistic, 
matter-dominated outflow reads (Bll) 



dlnr 



= (1-A*) 



-3/.)-(l-.^)|. 



r(i + A^) 



(7) 



' The radial component Br is suppressed as r ^ in the expanding plasma 
and can be neglected. 
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Here I^, is the specific intensity, is the photon frequency, 
fj, = cos 9, 9 is the angle relative to the radial direction, and 
r is the distance from the central source. The emission and 
absorption coefficients j^, and represent all the interactions 
of radiation with plasma as well as with radiation itself. All 
quantities (except r and F) are measured in the rest frame of 
the outflow. 

Equation (|7]l and the formulation of the ultra-relativistic 
transfer problem are discussed in detail in Bl 1. The problem 
simplifies because essentially all photons flow toward larger 
r, and hence only inner boundary conditions need to be spec- 
ified. The radiative transfer then takes the form of an initial- 
value problem. One can think of it as the evolution of with 
r or, equivalently, with the jet comoving time t. Equation (|7]i 
is similar to the kinetic equation for the evolution of the pho- 
ton distribution function (VP09), with additional terms due to 
the angular dependence. 

The inner boundary condition for is the blackbody radi- 
ation of a given temperature. We track the evolution of radi- 
ation with r consistently with the plasma evolution, which is 
described in Section lS^ In the case of a passively cooling jet 
(no heating), the transfer equation reproduces the adiabatic 
cooling of photons in the subphotospheric region. Heating 
changes the state of the plasma and greatly affects the source 
function Si, = jv/ Kd that governs the radiative transfer 

Our numerical solution of Equation d?) will use one ap- 
proximation to reduce the computational time: when calcu- 
lating ii, and we pretend that radiation is isotropic, using 
only the zeroth moment of the actual angular distribution of 
photons. In particular, the scattering rates are averaged over 
angles. The approximation is kno wn to be quite succes sful 
in classical transfer problems (e.g. IChandrasekharl IT960l) . It 
is also reasonable for GRB outflows if the bulk of scattering 
events are non-relativistic in the plasma rest frame (Bl 1). Its 
error becomes more significant at high photon energies, com- 
parable to mcC^ in the plasma frame. This approximation al- 
lows us to use the numerical tools for and k„ developed in 
VP09. A detailed discussion of the angle-averaged and 
and their numerical treatment is found in that work. 

Radiation that will be received by a distant observer is de- 
termined by at a large enough radius, where photons stream 
almost freely at all relevant v (note that this radius is much 
larger than iZ^). Using 1 — F^^ « 1, we define Doppler fac- 
tor as I? = F(l + The spectral distribution of the observed 
(isotropic equivalent) luminosity is given by the Lorentz 
transformation of d/i to the observer frame and integration 
over angles. 



L^{E) 



V{ti) dfi, (8) 



where E — Vhv is the observed photon energy and h is the 
Planck constant. We emphasize that the simple relation (O is 
valid only at large radii, after all transfer effects on have 
already been calculated according to Equation dTji. 

Finally, we note that Equation dTji was formally derived for 
steady outflows, however it is also applicable to strongly vari- 
able, non-uniform jets (see Bll and Appendix lAll. The basic 
reason for this is that radiation moves together with the ultra- 
relativistic plasma flow. The photon diffusion in the plasma 
rest frame is limited to scales 5 ^ ct ^ r/F, which are much 
smaller than r. When viewed in the fixed observer frame, the 
scale of diffusion in the radial direction is additionally com- 
pressed by the factor of F^^. The transfer occurs in a thin 



shell (or "pancake") of the jet as if it were part of a steady 
outflow. 

3.2. Kinetic equation for plasma 

The comoving-frame kinetic equation for electrons and 
positrons can be written as 

dn±{p) d 
dp 



dt 



1 d 

pn±{p) - - — (Dn±{p)) 



-J± 



K±n±{p) - - n±(p). 



(9) 



Here p = 1/7^ — 1 is the electron/positron momentum in 
units of JTieC n±{p) are the distribution functions of positrons 
and electrons, and t — r/cF is the proper time of the flow. 
When calculating p we take into account Coulomb heat- 
ing by protons (Equation |5]l and adiabatic cooling (Equa- 
tion [TO] below). The term p also includes the usual contri- 
butions from Compton scattering (in the Thomson regime), 
synchrotron emission and self-absorption, and Coulomb col- 
lisions between e^, as described in VP09. The diffusion coef- 
ficient D governs particle thermalization at low energies due 
to synchrotron processes and Coulomb collisions between 
(VP09). The source j± includes the injection of high-energy 
pairs by inelastic nuclear collisions as described in Section |2] 
The source and sink terms j± and k± also contain contri- 
butions from photon-photon pair production, pair annihila- 
tion, and Compton scattering (in the Klein-Nishina regime). 
Our code calculates all these processes using the exact cross- 
sections, as described in detail in VP09. Finally, the last term 
in Equation (|9]l takes into account the dilution of particle den- 
sities due to the two-dimensional (sideways) expansion of the 
coasting flow. 
The contribution of adiabatic cooling to p is given by 



Pad = 



3 t' 



(10) 



This is easy to verify by considering adiabatic expansion of 
J\fe monoenergetic electrons occupying a volume V. Their 
pressure — {N'c/V){p'^ /S'^jnieC^ determines the adia- 
batic cooling rate McmeC^d'y = — Pc dV . The coasting out- 
flow expands in two dimensions, V (x r"^ (x t^. Using 
dp = {'y/p)d'y, one gets Equation ( fTOl) . 

In our numerical simulations, the electron/positron kinetic 
equation (|9]l is discretized on a grid of particle momenta and 
solved simultaneously with the radiative transfer equation (|7]). 
The grid of dimensionless momentum p extends from 10^'' 
to 10*, with 25 points per decade. The photon-energy grid 
spans 15 decades from hv/rricC^ = 10~^^ to 10* (comoving 
frame), with approximately 13 gridpoints per decade. The 
lower boundary is set so low to properly simulate the effects 
of synchrotron self-absorption in mildly relativistic pairs in a 
weak magnetic field. The photon angular grid is uniform in 9 
and has 40 grid points. 

4. NON-MAGNETIZED OUTFLOWS 

We begin with the non-magnetized model that can be di- 
rectly compared with the results in BIO. The model has 
the following parameters: proton flow luminosity L = 
10^^ erg s^^, neutron flow luminosity Ln = 2 x 10^^ erg s^^, 
Lorentz factor of the proton flow (baryon loading) F = 600, 
Lorentz factor of the neutron flow F„ = 100, initial radius of 
the flow ro = 10^ cm. The starting radius of the simulations is 
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Figure 1. Momentum distribution function for electrons and positrons in 
tlie non-magnetized jet at different radii: r/Rn = 3, 10, 30, and 100 (solid, 
long-dashed, dashed and dot-dashed curves, respectively). The distribution is 
normalized so that it shows the optical depth rx , which is related to density 
ric = n+ + n- by Equation ^6). The dotted line indicates the slope of 
a nonthermal distribution that would be obtained if the injected did not 
produce a cascade and were Compton-cooled in the Thomson regime. 

determined by Equation dU, which gives = 10^^ cm. The 
comoving temperature of the blackbody radiation field at 
is found from the passively cooling outflow model at r < R^, 
which gives kTBB{Rn) ~ 0.5 keV. 

4.1. Electron/positron distribution 

At the radius R^^, where the dissipation begins, the injected 
high-energy pairs start to upscatter the thermal photons to en- 
ergies up to ~ lOOmoC^ in the rest-frame of the jet. The opti- 
cal depth for photon-photon pair production seen by the high- 
energy photons is initially small, as the target photons are con- 
fined to the blackbody component that is below the threshold 
for the reaction. However, this soon changes as more and 
more photons are upscattered from the thermal distribution. 
Then the jet becomes optically thick to pair production and 
a cascade develops. Soon the annihilation balance is estab- 
lished, in which the cooled pairs annihilate at the same 
rate as the new pairs are created by the cascade. 

The continual cascade quickly reaches a steady state, on 
a timescale shorter than the jet expansion timescale. Then 
a quasi-steady distribution is maintained, which gradu- 
ally changes with radius (Figure [T). The distribution has two 
parts, thermal and nonthermal (cf. BIO). The nonthermal 
tail consists of continually injected, fast-cooling pairs. If 
the injected generated no cascade and cooled in Thomson 
regime at all p, the nonthermal tail would have a power-law 
form ric(p) oc p^^. The actual distribution is affected by the 
Klein-Nishina effects in Compton cooling (which are signifi- 
cant atp > 50) and the large number of secondary pairs pro- 
duced in the cascade. 

The low-energy peak of the distribution maintains a 
Maxwellian shape due to Coulomb collisions between par- 
ticles. This thermalized population is continually heated by 
Coulomb collisions with protons and maintains the equilib- 
rium temperature that is determined by the balance between 
Coulomb heating and Compton cooling. The resulting elec- 
tron temperature is close to 15 keV and remarkably stable 
throughout the dissipation region (see Figure |4] below). Its 
weak radial dependence is approximately described by equa- 
tion (32) in BIO. 

Shortly after the dissipation starts, the scattering opacity of 



T 




E [MeV] 

Figure 2. Comparison of spectra obtained by kinetic (solid line) and Monte 
Carlo (dashed line) simulations for the non-magnetized jet. Model param- 
eters: proton flow luminosity L = 10^"^ erg s~^, neutron flow luminosity 
Ln = 2 X lO^'^ erg s~^, Lorentz factor of the proton flow F = 600, Lorentz 
factor of the neutron flow Fn = 100. The dotted straight lines correspond 
to photon indices 0.5 and —2.5. Note that the calculated spectra are not cor- 
rected for the cosmological redshift z. All photons should be redshifted by 
the factor {1 + z)~^ , which is 1/3 for a typical GRB. 

the flow becomes dominated by the created pairs. The 
Thomson optical depth is regulated to tt ~ 20 near R^ and 
then decreases as Rn/r (see Figure[T]l, in excellent agreement 
with BIO. This places the photospheric radius at i?^ ^ 20i?n- 
Since tt oc ricr cx r^^. Equation ^ shows that the 
Coulomb heating peaks near R^, and the fraction of the flow 
energy dissipated in one dynamical time at larger radii scales 
as r^^. The nonthermal energy injection has the same r^^ 
dependence. Thus, the ratio of thermal and nonthermal heat- 
ing rates remains nearly constant along the flow. This ratio is 
close to unity. 

Most of the kinetic energy of the flow is dissipated (and re- 
dissipated) by collisional heating. One should keep in mind 
though that the heat given to plasma and radiation in the sub- 
photospheric region is continually degraded by adiabatic cool- 
ing. As a result, the actual luminosity released at the photo- 
sphere is smaller by a factor of 1/2 compared with the total 
heat deposited in the flow (see section 4.6 in BIO). 

4.2. Radiation spectrum 

The spectrum that should be received by a distant observer 
from the collisionally heated jet is shown in Figure |2] For 
comparison, we also show the result of the Monte-Carlo sim- 
ulation in BIO. The two spectra are very similar 

The emitted spectrum peaks near 1 MeV. Radiation be- 
low the peak is made of blackbody photons advected from 
the central source and released at the photosphere. Note 
that the low-energy photon index a w 0.4 — 0.5 is signifi- 
cantly different from the Rayleigh- Jeans a = 1. This soften- 
ing is a result of the superposition of emissions from differ- 
ent parts of the photosphere having different Doppler factors, 
i.e. effectively a multi-temperature spectrum is observed be- 
low the peak (BIO). The spectrum between 1 and 20 MeV 
is mainly shaped by thermal Comptonization by Coulomb- 
heated pairs. The thermal pairs upscatter photons to energies 
up to £' ^ 2TkTa ^ 20 MeV in the observer frame, where 
kTc ^ Ih keV is the electron temperature in the comoving 
frame. 

Above 20 MeV, the observed spectrum is generated by in- 
verse Compton scattering by the nonthermal pairs. Note that 
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there is no dip between the thermal and nonthermal Comp- 
tonization components. Two reasons contribute to this re- 
markable smoothness of the spectral shape. First, the en- 
ergy budgets of thermal and nonthermal Comptonization are 
comparable, and the two components have comparable lumi- 
nosities at 10-20 MeV. Secondly, most of the photons above 
20 MeV undergo several scattering events on thermal pairs 
before escaping. They experience significant Compton down- 
scattering, which has an overall smoothing effect on the spec- 
trum. The only distinct feature in the predicted high-energy 
spectrum is a broad annihilation line on top of the smooth con- 
tinuum just below 1 GeV. The line is produced by the cooled 
pairs (photon energy E ^ 1 GeV in the observer frame 



Table 1 

Results of simulations. 



approximately corresponds to hiy 



in the jet frame). 



The observed spectrum at energies E > 1 GeV is affected 
by absorption due to photon-photon (7-7) pair production. 
The jet is optically thick to GeV photons close to where 
most of the heating occurs, and all high-energy upscattered 
photons are quickly absorbed, creating secondary pairs. 
Thus, the developing pair-photon cascade is initially in the 
saturated regime. 

The photon emissivity produced by the saturated cascade 
may be estimated analytically (Svensson 1987). An important 
parameter of the cascade is zt = (2/3)7oa;o, where 70 is the 
injection Lorentz factor of the primary nonthermal pairs and 
xq = hvo/meC^ is the typical energy of soft photons. The 
quantity zt determines whether the cascade takes place in the 
Thomson or Klein-Nishina regime (zt < 0.5 or > 0.5, re- 
spectively), the number of generations of secondary pairs etc. 
In our case a;o ~ 1/200 and 70 = 300, so zt ~ 1, i.e. the cas- 
cade starts in the Klein-Nishina regime at high energies and 
proceeds in the Thomson regime at lower energies. There are 
several generations of secondary pairs and photons, leading 
to a smooth overall spectrum. For the cascade with zt ~ 1 
the analytic solution pr edicts a flat emi ssivity of high-energy 
photons, i^ji, w const (ISvenssonlll987l) . The radiation inten- 
sity inside an opaque source is given by I^, = j^l ^u.^-y- Here 
the absorption coefficient ' and /3 is the photon 

index of the target radiation (typically photons of multi-MeV 
energy in the observer frame). Thus, the high-energy spec- 
trum oc is maintained inside the opaque jet. 

The escaping multi-GeV emission may be estimated by 
considering the evolution of 7-7 opacity with radius (BIO). 
The rate of nuclear collisions fades at large radii, but it still 
generates a significant cascade. Since the 7-7 optical 
depth decreases with r, the multi-GeV photons produced in 
the outer region have a chance to escape. The radius of 7- 
7 transparency at energy E scales approximately as R^^ oc 
E~^~^, and the power of the cascade scales as r~^. This 
leads to L^{E) oc E^^ (photon index j5 — 1)Q This rough es- 
timate suggests the change in photon index above a few GeV 
from /3 to /3 — 1. A similar steepening of the spectrum is 
observed in Figured Note that the high-energy radiation es- 
caping from the collisionally heated jet extends far beyond 
1 GeV, up to 100 GeV. 

Finally, note that our spectrum deviates from that of BIO 
more significantly in the GeV range. This is likely caused 
by our angle-averaged approximation for ji, and Ki, (Sec- 
tion l3.1b . Anisotropy effects are particularly significant at the 

^ Incidentally, the high-energy radiation inside the opac[uc source ij^ — 
Ji///^;y,77 also has the photon index (3 — 1. This coincidence isaresuU of the 
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"Magnetization defined in Equation (TT). 

''Peak energy of the EL e spectrum. 

"^Photon index in the 100-500 keV range. 

''Radiative efficiency e^ad = / L. 
Radius of the Thomson photosphere _R* relative to radius iJn = 10^^ cm 
where the dissipation starts. In the case of e-g, = 2, the jet remains optically 
thin to scattering throughout the collisionally heated region (i?* < Rn)- 

'Pair yield Y = M / 70 . where M is the secondary pair multiplicity and 
70 = 300 is the Lorenz factor of injected electrons. 

^Ratio of thermal and nonthermal heating rates at -R*/5 for cases with 
Eb ^ 10—2, and at the radius of maximum r-r for cases with eb ^ O-l- 

''Pair temperature at -R*/5 for £3 < 10- 2, and at maximum tt for eb > 
0.1. 

high-energy end of the spectrum that forms at large radii. At 
r ^ i?* the soft photon field is strongly collimated, and the 
scattering of nearly radial photons by relativistic pairs is pref- 
erentially backward. The backward-scattered photons see a 
higher 7-7 optical depth and less of them can get out. The 
anisotropy of 7-7 opacity is further enhanced by the very 
same collimation of the soft radiation field. These effects cre- 
ate particularly strong angular dependence of and k„ at 
high energies, which is missed by our angle-averaged approx- 
imation. 



5. MAGNETIZED OUTFLOWS 



Magnetization of the jet is described by the parameter 

eb = 



L 



L 



(11) 



particular scaling ju oc r and 
jet. 



in the collisionally heated 



where Ub = S^/Svr is the magnetic energy measured in the 
rest frame of the outflowing plasma. The quantity eb is the 
ratio of the Poynting flux L b (measured in the static frame) 
to the kinetic luminosity of the outflow L. We discuss below 
the GRB model with the same parameters as in Section|4]and 
investigate how the results change with increasing eb- We 
have calculated a set of models with sb — 10^"^, 10^^, 0.1, 
0.5 and 2. The corresponding ratios of (comoving) magnetic 
and radiation energy densities at the start of simulations are 
UBiRn)/U-y{Rn) = 0.006, 0.06, 0.6, 3 and 12. The results 
of our calculations are summarized in Table [T] and described 
in detail below. 

5.1. Electron/positron distribution 

Strong magnetic fields imply significant synchrotron cool- 
ing of the high-energy pairs injected by nuclear collisions, 
which can compete with Compton cooling. Since the syn- 
chrotron photons do not create secondary pairs, the magnetic 
field has a suppressing effect on the cascade. The multi- 
plicity of each subsequent pair generation in a saturated cas- 
cade is approximately proportional to Siad / (Crad +£b), where 
Eiad — Ly/L is the fraction of the jet energy carried by ra- 
diation. For magnetic fields approaching equipartition with 
radiation (i.e. ebL/L^ ^ 1), only two generations of sec- 
ondary pairs can make a significant contribution to the total 
pair multiplicity. 
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Figure 3. Momentum distribution of electrons and positrons at radius 
r/Rn = 3 in the magnetized jet. The jet parameters are the same as in 
Figures[T]and|2]except for magnetization. Different curves show models with 
different sb- sb = (solid), 10"^ (short-dashed), 0.01 (long-dashed), 0.1 
(dotted), 0.5 (dot-dashed) and 2 (triple-dot-dashed). The bump near p ~ 5 
is the result of synchrotron self-absorption, which becomes increasingly im- 
portant at high eb . The dotted line indicates the slope of a nonthermal dis- 
tribution that would be obtained if the injected lost energy only by emit- 
ting synchrotron radiation, with no inverse-Compton cascade and no self- 
absoiption of synchrotron radiation. 

As a result, the total pair yield decreases from Y « 0.2 to 
about 0.005 as eb increases from to 2 (Table [T). The pair 
yield is defined as y = A4/jo, where Ai is the multiplic- 
ity of secondary pairs. In the model with £b = 2, the pair 
yield is not far from the minimum value Y = I/70, which 
corresponds to = 1 (no secondary pairs). The strong field 
completely suppresses the cascade. 

Figure Q\ shows the distribution function at r = 3i?„. 
The overall shape of the distribution is similar to that found 
in non-magnetized jets. It has thermal and nonthermal parts. 
With increasing eb, the normalization of the distribution sig- 
nificantly decreases — the number of is reduced as the 
cascade is suppressed. At eb ^ 0.1 the synchrotron cooling 
becomes dominant and the e distribution at high energies 
is described by the power-law nc{p) oc p^^. A new feature 
appears when the magnetization is high — the bump in the 
electron spectrum at p < 10. This bump is a result of syn- 
chrotron self-absorption that tends to "thermalize" the high- 
energy pairs before they join the thermal population at low 
energies and thermalize via Coulomb collisions. A similar 
effect is seen in the spectral simulations of accreting black 
holes in X-ray binaries and active galaxies (iPoutanen & Vumil 
12009; Veledina et al. 2011). 

The reduction in pair production implies a smaller Thom- 
son optical depth (tt oc Y^^^, see equation 23 in BIO). As 
a consequence, the photospheric radius i?* decreases with in- 
creasing Eb (Table Another significant implication is the 
reduction in the thermal heating rate (jth, which is propor- 
tional to Tt (see Equation|5]l. Since the nonthermal injection 
rate (jnth remains unchanged, the ratio (jth/Onth is reduced 
in magnetized jets. 

The temperature of thermalized pairs remains remarkably 
stable as eb is increased (see Figure |4). This fact may be 
understood by noticing that the thermal heating rate per par- 
ticle does not decrease with increasing magnetization (the 
volume heating does). Note also that the thermal pairs are 
unable to cool via synchrotron emission because of strong 
self-absorption. As a result their equilibrium temperature is 
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Figure 4. Temperature of the thermalized population as a function of 
radius, Tc{r). The jet parameters are the same as in Figure|2] except for mag- 
netization. Different curves show models with different eb - £b = (solid), 
IQ-^ (short-dashed), 0.01 (long-dashed), 0.1 (dotted), 0.5 (dot-dashed) and 
2 (triple-dot-dashed). At r/Rn < 10 the temperature is set by the balance 
between Coulomb heating and Compton cooling; in this region Tc{r) is well 
described by equations (32) and (37) in BIO. At larger radii, the contribu- 
tion of adiabatic cooling to the thermal balance becomes non-negligible. As 
a result, Tc begins to slowly decrease. 




E [MeV] 

Figure 5. Radiation spectrum emitted by the magnetized, collisionally 
heated jet. The jet parameters are the same as in Figure|2] except for magne- 
tization. The solid, short-dashed, long-dashed, dotted, dot-dashed and triple- 
dot-dashed curves correspond to magnetizations eb = 0, 10^'', 0.01, 0.1, 
0.5 and 2, respectively. The straight dotted line shows a power-law spectrum 
with a = —1.2. The spectrum L^; is multiplied by photon energy E to make 
the differences between the models more visible in the figure. 

weakly affected by the magnetic field. 

5.2. Radiation spectrum 

The spectrum that should be observed from the magnetized 
jet is shown in Figure|5] Magnetization significantly changes 
the radiative properties of the jet. It suppresses the pair cas- 
cade, reduces the photospheric radius, and changes the dis- 
tribution function. It also implies a new emission component 
— synchrotron emission from the nonthermal pairs. 

Our detailed calculations confirm that the synchrotron radi- 
ation from strongly magnetized collisionally heated jets peaks 
at energies comparable to 1 MeV in the observer frame, as 
estimated in previous works (Koers & Giannios 2007; BIO). 
The position of the synchrotron peak varies with magneti- 
1/2 

zation as £3' . Its amplitude increases linearly with £3 for 
weak magnetizations (£3 < 0.01) and approaches a constant 
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value when synchrotron emission becomes the dominant en- 
ergy loss mechanism for high-energy pairs. The synchrotron 
peak remains practically buried under the Comptonized ther- 
mal spectrum in all models shown in Figure |5] Thus, the 
emerging spectrum preserves the Band-type shape with the 
MeV peak. 

The presence of synchrotron radiation significantly affects 
the low-energy slope of the Band peak — it makes the spec- 
trum softer, i.e. a is significantly reduced compared with the 
non-magnetized model. The obtained spectral shape below 
1 MeV may be described as follows. The mixture of syn- 
chrotron and Comptonized thermal photons create a nearly 
power-law spectrum between 100 and 500 keV (with a slope 
a that depends on eb)- Below this energy range, the spec- 
trum is dominated by synchrotron radiation that has a smaller, 
softer slope. The resulting curvature of the spectrum may be 
described by observers as a soft excess above the power law. 
At energies < 10 keV the spectral slope changes again, 
as s ynchro tron self-absorption becomes important (see Sec- 
tion l5T2] l. 



5.2.1. High-energy emission 

Three mechanisms contribute to emission at energies E > 
1 MeV: thermal Comptonization by pairs with temperature 
~ 15 keV in the jet frame, nonthermal Comptonization by the 
cascade, and synchrotron emission (which extends to tens 
of MeV when e-q is large). 

Thermal Comptonization dominates at i? > 1 MeV, near 
the spectral peak. The magnetized jets have approximately 
the same electron temperature Tc as non-magnetized jets, 
and their optical depth tt is smaller As a result, the Kom- 
paneets' parameter y — ATThTc/nicC^ is reduced with in- 
creasing eb- This leads to a steeper slope of the thermally 
Comptonized spectrum above the peak. Similar to the non- 
magnetized model, the thermally Comptonized power law de- 
clines at > 2rfcTo - 20 MeV. In models with high £b, 
synchrotron emission makes a comparable or even dominant 
contribution at these energies. In the model with £b = 2, syn- 
chrotron emission contributes significantly to the spectrum up 
to 50 MeV. 

The luminosity above 100 MeV is produced only by in- 
verse Compton scattering by the nonthermal particles. This 
luminosity is inevitably reduced with increasing eb, as part of 
the energy of injected pairs is lost to synchrotron emission at 
lower energies. The nonthermal inverse Compton component 
becomes weaker and harder with increasing £b, and its slope 
approaches Ly{E) oc E^^^^. This slope is the signature of 
inverse-Compton emission with suppressed pair cascade. 

Overall, the suppression of the pair cascade by synchrotron 
cooling destroys the simple power-law shape of the high- 
energy spectrum. Instead, a distinct hard component (non- 
thermal inverse Compton) appears above 50-100 MeV. 



5.2.2. Low-energy emission 

The low-energy end of the predicted spectrum is dominated 
by synchrotron emission, even when £b is small (Figure |5]l. 
The spectrum at energies E < 1-10 keV is affected by self- 
absorption. It can be derived analytically as follows. 

In the rest frame of the plasma, the angle-averaged syn- 
chrotron emissivity and the absorption coefficient are given 



by (see e.g. lGhisellini & Svensson|[r99Ih 



js{v,p) Jlcip) dp, 



dp 



(12) 
dp. (13) 



Here all quantities are measured in the plasma rest frame; 
nc{p) — n+{p) + n-{p) is the distribution function of 
pairs, and js{i^,p) is the angle-averaged synchrotron emissiv- 
ity per electron. For analytical estimates we will use the delta- 
function approximation for the emissivity 

4:TTjs{i^,p) = ^ccttUb p^ 5(y - j'^i^b) , (14) 

where i'b = eB /2TimcC is the Larmor frequency. The syn- 
chrotron emission is produced by relativistic particles with 
7 w p. Then Equations (fT2b-(fT4l) give 



is{v) = — hvBpnc{p), 



Ks(j/) 2 + 5 



nici'BP^, 



(15) 



(16) 



where a/ = e^/hc — 1/137, 6 — —d\nnc{p)/d\np is the 
local slope of the distribution function, and 

P«7-( — ) ■ (17) 

The distribution function ndp) in Equation (flSl l can be de- 
termined by assuming a quasi-steady flow of particles in 
the momentum space and writing 



pucip) 



_ - inj 



M{p), 



(18) 



where n™^ is the rate of particle injection at the highest en- 
ergy 7o ~ 300 (Equation O and M{p) is the multiplicity of 
secondary pairs created with momenta above p. The syn- 
chrotron energy losses for particles emitting in the optically 
thin regime are given by 



2 £b 4 . 

7moC ; = -c(JtUbP 

£rad + Eb O 



Using Equations (flSb and ( fT9] l (with p 
emissivity (fTSl l becomes 



(19) 

7) the synchrotron 



87ri/B(erad + £6)7' 



(20) 



Let's now evaluate the range of Lorentz factors 7 > 7s for 
particles that emit synchrotron radiation in the optically thin 
regime, as a function of radius r. The synchrotron photo- 
sphere can be found from the approximate condition 



1. 



(21) 



Using Equations (fTSl l and ( l20b . together with the relations (O 
and (fTTI ). we find from Equation 



3(2 + S)eb 



07 / , ^ ^n^TnM{-rs)- (22) 
2^7r(£rad + sb) r„ 
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With [2(2 + (5)] 1/*^ ft 
(me/mp)i/3 



23/4 1/12 j.l/3 



2^/-^, this equation gives 

(ffrad + eB)<JTLVl 



-,1/6 









1/12 

(23) 



where Tc = jm^c^ ~ 2.82 x 10^^'^cmis the classical elec- 
tron radius. This gives 7s ~ 10 for the parameters adopted in 
our models. Note that 7s does not depend on r and weakly 
depends on the parameters of the jet. 

Synchrotron radiation at a given radius r is self-absorbed at 
energies E < where 



E,{T)^Y^lhvB- — ^^ 
rricC r 



IebL 



1/2 



(24) 



The entire heating region r > i?n is transparent to syn- 
chrotron absorption for photons above Es{Rn)- For the jet 
models calculated in this paper Es{Rn) ^ 1-10 keV. Ob- 
served radiation at lower energies E comes mainly from the 
corresponding synchrotron photosphere Rs{E) that can be 
found by expressing r from Equation ( l24l i. 

Neglecting the self-absorbed radiation from r < Rs{E),we 
estimate the observed spectral luminosity Ls{E) by integrat- 
ing the synchrotron emissivity over the optically thin region 
r > RsiE), 



Ls{E) 



Anr dr. 



(25) 



Here we approximated hiy w E/T and used the fact that the 
angle-integrated emissivity Anjs is the same in the static and 
comoving frames within a factor C ~ iQ Using the slope 
of distribution (5 « 2, one can show that js{v)r^ scales 
with radius as r^'^^+'^^Z^ = 7--3/2 ^jj^j jj-g integral peaks near 
Rs{E). Substituting js{v) from Equation (|20] |. using the re- 
lations 7^ ^ E/ThvB, r^B^r^ = 2eBi/c and Equation ©, 
we obtain 



L 



rricC^ 



r-^7,^eB, E<Es{Rn), (26) 



where 7s is given by Equation ( |23] ). The emission at any pho- 
ton energy E < Es{Rn) peaks near Rs{E) and is produced by 
particles with the same 7s ~ 10. Note that Ls{E) = const 
(flat spectrum), which corresponds to photon index a = — 1. 
This fact is a consequence of B oc r^^ and nc{p) oc r^^. For 
a similar reason , a flat spectrum was der ived for opaque radio 
jets in AGN (.Blandford & K6nig]|[T979t) . Equation ^ is in 
approximate agreement with our numerical results (Figure 0. 

5.3. Radiative efficiency 

The radiative efficiency of the jet ei ad is defined as the ratio 
of the photon luminosity i-y to the kinetic luminosity of the 
plasma outflow L. Note that L{r) « const for the matter- 
dominated jet and evolves with radius. This evolution is 
shown in Figure|6] The final efficiency is given by the asymp- 
totic value of Eiad = L^/L at large radii. 

The photon luminosity prior to the onset of dissipation 
(r = i?n) is found from the passively cooling jet model. One 



^ This factor is given by 



C ■■ 



1 /D 

2 7-1 Vr 



1 "+ 
2 J ^(1 + ^)"'^^' 



5+1 




0.0 



passive 



10 



10^ 



r/R„ 



Figure 6. Fraction of the flow energy carried by radiation (radiative effi- 
ciency) as a function of distance from the central source, for different mag- 
netizations: eb = (solid line), 0.01 (long-dashed), 0.1 (dotted), 0.5 (dot- 
dashed) and 2 (triple-dot-dashed). For compaiison, the radiative efficiency 
for a passively cooling flow is also plotted (lower solid line). 

can see from Figure |6]that is greatly increased by the in- 
tense collisional heating at r > i?n- The net energy given to 
radiation by collisional heating outside a given radius is ob- 
tained by integrating Equations (|3]l and ^ over volume; the 
result is proportional to r^^. Thus, heating peaks near Rn 
and continues with a smaller rate at larger radii. On the other 
hand, is reduced by adiabatic cooling, in particular in the 
opaque zone r < Ri,. The competition between collisional 
heating and adiabatic cooling shapes L^ir). In weakly mag- 
netized jets, the photospheric radius is large, i?^ ^ 20i?n, 
and adiabatic cooling is more efficient. It begins to win over 
collisional heating at r ^ 5i?n and somewhat reduces L^. 
In strongly magnetized jets, the photospheric radius i?* is 
smaller because the magnetic field suppresses the production 
of pairs (Section lSTt . In this case, adiabatic cooling is less 
efficient. 

In all cases shown in Figure |6] the final Crad is close to 
50%. We conclude that the radiative efficiency of collisionally 
heated jets remains high for the entire range of Eb considered 
in this paper 

6. DISCUSSION 

Only two radiative processes are in principle capable of pro- 
ducing the GRB spectrum with an extended high-energy tail: 
synchrotron emission and Comptonization. These are two ba- 
sic modes of energy transfer from a heated plasma to radia- 
tion, and their theory was developed long ago. In particular, 
the formation of photospheric spectra through Comptoniza- 
tion in X-ray sources was investigated in the 1970s, with var- 
ious applications from accretion disks to the expanding Uni- 
verse (see e.g. Pozdnyakov et al. 1983). Comptonization in 
GR Bs h as been discussed for about two decades (see Sec- 
tion l6.2l below). Besides the fact that subphotospheric heating 
creates a nonthermal-looking radiation spectrum, one would 
like to know what the heating process is. The answer to this 
question can help disentangle GRB physics — the nature and 
composition of the jet and the central engineQ Progress in this 
direction was hampered for many years by the complexity of 
collisionless processes that were usually invoked in GRB pro- 

* The situation may be compared with X-ray binaries, where Comptoniza- 
tion is thought to occur in a corona of the accretion disk or the inner hot flow, 
and the main puzzle is how the plasma is heated. 
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duction. 

GRB spectra can be shaped by three possible heating mech- 
anisms: (1) collisionless shocks, (2) magnetic dissipation, and 
(3) collisional dissipation. Least understood is magnetic dissi- 
pation, and it is usually modeled by introducing phenomeno- 
logical parameters. A similar phenomenological approach 
was used for internal shocks in the jet until recent numeri- 
cal simulations began to provide insights into shock physics. 
Collisional dissipation is the most straightforward mechanism 
of these three. It can be calculated exactly, from first princi- 
ples. BIO recently showed that collisional heating possesses 
the key features of a successful GRB mechanism: efficient 
electron heating, high radiative efficiency and, most impor- 
tantly, it generates the observed Band-type spectrum. 

In this paper, we calculated the emission from coUisionally 
heated jets by solving the time-dependent, coupled kinetic 
equations for the particle and photon distributions inside the 
jet. The advantage of our numerical code is the accurate mod- 
eling of all relevant kinetic and radiative processes including 
Coulomb collisions and synchrotron self-absorption. It al- 
lowed us to systematically explore the effects of jet magneti- 
zation on the emerging spectrum. We calculated the emerging 
radiation using the equation of radiative transfer (B 11 ) instead 
of the Monte-Carlo method used by BIO. The only disadvan- 
tage of our code is that it uses angle-averaged coefficients of 
emission and absorption in the transfer equation. Compari- 
son with exact transfer models of BIO and Bll shows that 
this approximation is reasonable, in particular at photon en- 
ergies below GeV. Note also that this is the first implementa- 
tion of radiative transfer calculations in a kinetic GRB code. 
The previously developed kinetic codes ass umed isotropy of 
radiation in the rest-frame of the plasma ( Pe'er & WaxmanI 
120051: VP09), which is invalid in the main region of interest 

TT < 10. 

6. 1 . Spectra from non-magnetized and magnetized jets 

First, we calculated the emission from the fiducial baryonic 
jet model of BIO with zero magnetic field (eb = 0). Our re- 
sults agree with BIO: the spectrum peaks at 1 MeV, has the 
low-energy photon index a w 0.4 and the high-energy photon 
index (3 w —2.5. We conclude that the Band-type spectrum is 
a robust prediction of the collisional-heating model. Similar 
to BIO, we stress the ab initio character of our calculations: 
the dissipation process and generated radiation are derived 
from first principles, without introducing any phenomenolog- 
ical parameters for the heating mechanism. 

The collisional heating peaks below the photosphere, where 
the 7-7 opacity is large for multi-GeV photons emitted in 
the cascade. The nuclear collisions also occur at much 
larger radii (although with a smaller rate) where the 7-7 op- 
tical depth is reduced and the high-energy photons can es- 
cape. As a result, instead of a cutoff at a few GeV, the pre- 
dicted spectrum exhibits a moderate steepening at ^ 5 GeV 
and extends up to ^ 100 GeV (Figure |2]i. A similar behav- 
ior is seen in the Monte-Carlo results of BIO, with somewhat 
smaller normalization of the multi-GeV emission. This dif- 
ference must be caused by the angle-averaged approximation 
for ji, and that was adopted in our calculations of radiative 
transfer. As discussed above, the accuracy of this approxima- 
tion is reduced at photon energies E > 1 GeV. 

Then we calculated the emission from magnetized jets. We 
considered the range of magnetization parameters < Eb < 
2, so the magnetic energy was at most comparable to the en- 
ergy of the baryonic component. We did not consider mag- 



netically dominated jets with eb ^ 1 for two reasons. First, 
the strong magnetic field can change the jet dynamics. Then 
the coasting approximation (which is reasonable for a matter- 
dominated jet) would be invalid. Second, the strongly mag- 
netized case is computationally more difficult. The generated 
spectra in this regime are left for a future study. 

Our calculations show that magnetization has a small effect 
on the position of the spectral peak -Epoak — it is only slightly 
shifted to lower energies. The spectrum still has the Band- 
type shape around the peak, with a and (3 depending on e-q 
(Figure|5]and Table[T]i. Magnetization significantly affects the 
shape of the spectrum at high {E 3> i?pcak) and low {E <C 
-Epcak) energies. 

The high-energy emission is reduced with increasing e-q 
(Figure |5]l, because the nonthermal generated by nuclear 
collisions emit less via inverse Compton scattering and more 
at low energies via the synchrotron mechanism. This change 
is significant if eb S> 10^^. Synchrotron cooling of the in- 
jected high-energy has a throttling effect on the pair cas- 
cades. As a result, magnetized jets have smaller densities of 
pairs and smaller photospheric radii i?*. 

The low-energy emission is increased with increasing eb- 
An important effect of magnetization is the softening of the 
spectral slope below i?pcak- This is the result of synchrotron 
emission from nonthermal pairs. At energies below ^ 30keV 
a clear soft excess is predicted, even for small eb ^ 10"'^. 
Our model also predicts that the spectrum below a few keV 
should have the photon index —1. This spectrum extends 
down to the optical band. 

Note that photon spectra in this paper are not corrected for 
the cosmological redshift z. The actually observed energy of 
all photons is smaller by the factor of (1 + z)"^. In particular, 
the soft excess is redshifted to ^ 10 keV for a typical z of 
GRBs. 

6.2. Comparison with other models of subphotopsheric 
Comptonization 

iThompsonl (119941) proposed that GRB photons are Comp- 
tonized in the subphotospheric region by turbulent bulk mo- 
tions of the plasma in the jet. His model assumed that the jet 
is dominated by magnetic field and pictured Alfvenic turbu- 
lence generated by some instabilities (e.g. by reconnection). 
Building a detailed theory of turbule nce from firs t principles 
would be a formidable task, and Thompson (1994) made esti- 
mates assuming that the plasma motions are limited by radia- 
tion drag. Then the bulk Comptonization occurs in the unsat- 
urated regime and controlled by the efficiency of turbulence 
generation against the drag. The predicted GRB spectrum is 
sensitive to this unknown efficiency, which may vary with ra- 
dius. The mechanism is very different from collisional heat- 
ing. Direct comparison of the model predictions is difficult, 
because Thomp son ( 1994 ) used the theory of Comptonization 
in static sources, which is not valid for GRBs (see discussion 
in Section[3]). Like in the expanding Universe, the number of 
scatterings in a GRB jet scales as ~ tt (not as t^) and is con- 
trolled by the decrease in density, not the escape of photons 
from the plasma. The expected spectrum can only be found 
by solving radiative transfer in the expanding jet, and we em- 
phasize the importance of accurate transfer calculations for 
GRB models 

^ When heating occurs in a static source, the emerging spectmm can have 
two peaks, resembling the Comptonized spectra of accretion-disk corona. In 
contrast, the spectrum emitted by the GRB jet has one peak. 
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iMeszaros & ReesI (l2000lj) and iRees & MeszarosI (l2005h 
discussed subphotospheric Comptonization in general terms 
and highlighted its importance for GRBs. They argued that 
Comptonization helps explain observations, although they did 
not consider any co ncrete mod el. More detailed calculations 
were performed by PeLetet al. (2006). They still did not fo- 
cus on any concrete physical scenario and experimented with 
various phenomenological models, using free parameters to 
describe the heating process: the amount of dissipated energy 
and its fraction that is given to electrons, the shape of the elec- 
tron distribution, and the optical depth at which the dissipation 
occurs. Various combinations of these parameters led to vari- 
ous spectra. In contrast, we investigated the concrete physical 
model of coUisional dissipation, which predicts the thermal 
and nonthermal heating rates, both scaling as r~^. It is there- 
fore hard to compare their r esults with ours. We can only 
compare the technical tools. Pe'er et al.' ('2006) used the ki- 
netic code developed by Pe'er & Waxman (2005). In many re- 
spects, their code is similar to ours. Their code assumes, how- 
ever, that radiation is isotropic in the comoving frame of the 
jet. We gave up this approximation because it is violated even 
below the photosphere (see Section O and instead solved the 
radiative transfer equation. Note also that Pe'er et al. (2006) 
argued that GRB jets never become dominated by pairs. 
In contrast, we find that jets with magnetization eb ^ 1 have 
a large pair loading factor rio/n 20 — 40, which changes 
the photospheric radius by the factor of 20^0. 

Giannios (2006) investigated Comptonization in jets that 
carry alternating magnetic fields and are heated by magnetic 
reconnection (D renkhahn & Spruit 200 2). This model has 
two phenomenological parameters — the reconnection rate 
and the fraction of dissipated energy that is deposited into 
electrons. The effect of magnetic dissipation on the electron 
distribution is unknown. Giannios (2006) assumed a ther- 
mal electron distribution and studied the radiation produced 
by this model. He used a Monte-Carlo code and an iterative 
technique to simulate the radiative tran sfer in the jetQ We note 
that the jet model in iGiannios I (l2006h is significantly differ- 
ent from BIO and our work. In particular, the jet continues to 
accelerate as F cx r^/'^ in the heated subphotospheric region 
while we focus on matter-dominated jets with F « const. 
The electron temperature scales with radius approximately as 
Tc oc r^/'^ while in our model we find that To varies only by a 
fact or of ^ 2 f or three decades in radius (Figure 2) • Note also 
that IGiannios I (f2006) assumed pure thermal heating, leading 
to negligible pair creation. In contrast, our work and BIO pre- 
dict a nonthermal cascade, which has significant effects on 
the observed spectrum at high energies (above 10 MeV) and 
at low energies where synchrotron emission from the cascade 
b ecomes dominant. 

iLazzati & Begelman I (120101) considered a simplified model 
for Comptonization following an impulsive heating event be- 
low the photosphere. In their model, all electrons in a slab 
of optical depth tt 2 are suddenly given a Lorentz fac- 
tor 7 and then cooled by blackbody radiation. This setup can 
produce a variety of photon spectra depending on photon-to- 
electron ratio n^/rie, tt, and 7. It has, however, a problem. 
The sudden heating of the slab (faster than Compton cool- 
ing) can hardly describe any physical dissipation mechanism, 

* The numerical methods used by Giannios 1 420061) and BIO a re similar, 
except that the iterations of the electron temperature in lGiaimios I (S)06) are 
simplified: the profile Tc(r) is searched in the power-law form, with two 
parameters to iterate — the normalization and the slope of the power law. 



as different parts of the slab are causally disconnected. The 
causal contact (light-crossing) time for a region with tt ~ 2 
is orders of magnitude longer than Compton cooling time. 

6.3. Comparison with observed GRB spectra 

The main part of observed GRB radiation is emitted with 
a Band-type spectrum — a smoothly broken power law that 
peaks near 0.3 MeV. The spectral slope s below and above the 
peak (a and /3) vary (iPreece et al.ll200()h . The most frequently 
measured low-energy slope a is close to —1, but in some cases 
it reaches 1 (Rayleigh-Jeans slope of the Planck spectrum). 
The most frequently measured high-energy slope /3 is near 
—2.5; it also significantly varies from burst to burst, roughly 
between —2 and —3. 

The spectra predicted by our model appear to agree with ob- 
servations. The predicted spectral peak is close to MeV/(l + 
z). The average observed values of a ^ — 1 and /3 —2.5 are 
consistent with the collisionally heated jet with magnetization 
Eb ^ 10^^ — 10^^. A lower magnetization eb < 10^^ leads 
to harder a up to 0.4. Larger slopes are not expected from the 
matter-dominated jets. The slope a = 1 is predicted only in 
the radiation-dominated regime - in this case the jet emits a 
Planckian spectrum (Bl 1). 

It should be noted that our calculations give practically in- 
stantaneous observed spectra. The main part of the prompt 
emission (photon energies from soft X-rays to a few GeV) is 
dominated by the photosphere in our model. The emission 
is produced by a sequence of short-lived (in observer time) 
independent emitters passing through i?^, — the sequence of 
shells of thickness 6r < Ri./T'^, which corresponds to dura- 
tion (5fobs ^ 10^'* s for typical parameters of GRB jets. The 
spectra obtained by modem detectors do not provide this high 
temporal resolution, because of insufficient photon statistics. 
The observed spectra must be a mixture of emissions from 
different shells, possibly with significant variations in L and 
^'pcak- The superposition of different instantaneous spectra 
generally tends to reduce the observed slope a and can give 
a ^ —1 even for non-magnetized jets (R. Mochkovitch, in 
preparation). 

An interesting feature of the magnetized model of colli- 
sionally heated jet is the soft excess that appears at ener- 
gies ~ 30(1 + z)^^ keV. It may provide an explanation 
for the X-ray excesses that have been observed in several 
bursts (Preece et al. 1996i), includ ing the recent example of 
GRB 090902B ( Abdo et aLl l2009al) . 

Prompt optical emission is expected from collisionally 
heated jets, with a slope L^{E) = const (photon index —1). 
This emission peaks at the self-absorption radius Rs{E) ^ 
Ri, and should be observed with a time lag ^ Rs/cT^ ~ 1 s 
with respect to the MeV radiation that is released at r 
i?*. The prompt optical emission has been detected in some 
GRBs, but the data are sparse, as prompt optical observations 
are difficult. In some cases (e.g. GRB 0803 19B) the optical 
luminosity appears to exceed the predictions of our model, 
suggesting an additional source of emission (e.g. shocks or 
neutron decay, see BIO). 

The high-energy index (3 predicted by the collisional model 
can easily be steeper than —2.5, because of a large e-g, (Fig- 
ure |5]l or less efficient heating (BIO). The model can also 
accommodate very hard slopes (3 > —2 that are observed in 
some GRBs (see e.g. Figure 7 in BIO). 

Our model predicts significant multi-GeV emission, espe- 
cially when £b is small. It implies that the prompt GRB emis- 
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sion can contribute to the radiation observed by LAT instru- 
ment of the Fermi telescope. The prompt high-energy photons 
should mix together with radiation from the other, long-lived, 
high-energy source that was identified by LAT in a fraction of 
GRBs and whose origin is debated. 

Overall, our results support the view that the main peak of 
GRB spectra is dominated by photospheric radiation of the 
ultra-relativistic jet from the central engine. Note that pho- 
tospheric emission is often erroneously pictured as a Planck- 
ian component in GRB spectra. As shown in Bll, the pho- 
tospheric spectrum has a Planckian shape only if the jet en- 
ergy is strongly dominated by radiation at i?*. This radiation- 
dominated regime requires an extremely high Lorentz fac- 
tor, which corresponds to a very baryon-poor jet. Baryon- 



dominated jets cannot produce Planckian spectra, even if the 
jet experiences no dissipation and only emits adiabatically 
cooled thermal photons from the photosphere (Figure 1 in 
BIO). In reality, the jet should experience strong coUisional 
dissipation, which has a huge impact on its photospheric ra- 
diation and makes it even more different from Planckian. 
The only feature shared with Planckian radiation is the single 
spectral peak at i?pcak ^ 1 MeV. The photospheric spectrum 
of the collisionally heated jet is much broader than Planckian 
and fits the observed GRBs. 
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APPENDIX 
RADIATIVE TRANSFER 



Equation (|7]) is valid for relativistic outflows in the limit F cxd. For finite F ^ 1, the equation accurately describes the 
transfer in a steady jet. Its accuracy is reduced if the jet is variable on scales 5r < r/V"^, as discussed in Bll. Here we give 
a more formal discussion of transfer in variable jets and the applicability of the steady equation (|7]). The general equation of 
radiative transfer in a spherically symmetric outflow with any Lorentz factor F(i, r) can be written as (iMihalaslll980h 



CI, 



d 



9 In J/ 



where C is the differential operator 



^^^(' + ^'^) J^+^('^ + ^)|:■ 



(Al) 



(A2) 



Here /? = (1 — F~^)^/^ is the flow velocity in units of c, t and r are time and radial coordinate measured in the static frame, 
whereas the cosine fi of the photon angle relative to the radial direction, the specific intensity I,^, the emission and absorption 
coefficients ji, and k^, are all measured in the local comoving frame. 
Let us define a coordinate s by the following implicit relation 



ct 



(A3) 



For a matter-dominated (coasting) flow, F(s) w const and s can serve as a Lagrangian coordinate that labels different shells of 
the flow, as s = const along the worldlines of the shells. Changing variables t ^ r, s, we rewrite the operator (IA2) as 



C = T{p + /3) 



d_ 

dr 
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1 



dlnF 
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Here the spatial derivative d/dr is taken along the worldline (i.e. at s = const). Note that the flow worldlines do not cross (i.e. 
there are no caustics) as long as the bracket in Equation (IA4t is non-zero^ 
If the flow varies on scales 5r ^ r/F^, the two parts of the operator (IA4) are comparable. The problem simplifies (becomes 

formally equivalent to the steady problem) if the flow varies on scales much larger than r/F^. Then the operator C is dominated 
by the first term that is proportional to d/dr. The left-hand side of Equation ( lAU simplifies to CI, = r(/i + /?) dly/dr. 
Similarly, the right-hand side of Equation (lAlt simplifies and the transfer equation becomes, 

d r, . 9 ^ 1 r , . ^ N . 9s dl, d 



1 



r-2 5 In ? 



dhii 



dfi 



[{^~^i^){l + ^3^i)I, 



(A5) 



which becomes Equation ^ for outflows wi th /3 « 1. Formally, this is achieved simply by replacing d/dr\t by the derivative 
along the flow worldline d/dr\s. Equation (IA5b can also be written using the flow proper time as the independent variable, 

t = r/cPV. 
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